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The action of weak, streamwise vortices on a plane, incompressible, steady 
mixing layer is examined in the large-Reynolds-number limit: The outer, 
inviscid region is bounded by a vortex sheet to which the viscous region is 
confined. It is shown that the local linear analysis becomes invalid at 
streamwise distances 0(e _1 ), where e « 1 is the crossflow amplitude, and a 
new nonlinear analysis is constructed for this region. Numerical solutions of 
the nonlinear problem show that the vortex sheet undergoes an 0(1) change in 
position and that the solution is ultimately terminated by the appearance of a 
singularity. The corresponding viscous layer shows downstream thickening, but 
appears to remain well-behaved up to the singular location. 

1. Introduction 

Mixing layer flows being canonical and of technological importance have 
been studied extensively over the years. In this paper, some new, three- 
dimensional and nonlinear features are examined theoretically. We consider a 
steady, incompressible free shear layer--two nominally uniform streams at 
different speeds--and study the effect of small nonuniformities imposed on one 
of these streams. The analysis is restricted to the asymptotically large- 



Reynolds -number case for which the usual boundary layer approximation applies. 
There is an outer, inviscid region with an infinitesimally thin vortex sheet 
forming the interface between the two streams and an inner, viscous layer that 
resolves the velocity discontinuity across the vortex sheet. (See figure 1.) 

We suppose, for definiteness, that the mixing layer forms downstream of an 
infinitesimally thin splitter plate. The upstream base flow consists of a 
uniform parallel stream and a second stream at rest in our reference frame. 

The resulting vortex sheet is then a plane surface in the vicinity of the 
trailing edge, and the relevant viscous problem admits a similarity 
solution 1 ' 2 there. (For simplicity, we neglect the upstream boundary layer on 
the surface of the plate.) We suppose that a small 0(e) streamwise vorticity 
field is imposed on this flow somewhere upstream of the trailing edge, where 
e « 1 measures crossflow amplitude, and seek a series solution by expanding 
in powers of e. The 0(e) crossflow is convected without change, and the 
vortex sheet remains planar to that order. However, the 0(e 2 ) solution 
increases linearly with the streamwise coordinate x as x becomes large, so 
that the series is no longer valid when x is 0(l/e). A new, nonlinear 
solution therefore has to be found for this region in order to obtain a 
uniformly valid solution to the problem. A similar situation was found in a 
study of flat-plate boundary layers subjected to weak nonuniformities in the 

external stream. 3,4 

As in Goldstein and Leib, 4 the nonlinear flow is governed by the two- 
dimensional, time-dependent vorticity equation with an appropriately scaled 
streamwise coordinate playing the role of the time. However, the situation is 
now considerably more complex because the vortex sheet is displaced by an 0(1) 
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amount in the nonlinear region, and the vorticity equation therefore has to be 
solved in a region of unknown shape. The problem is also similar to those in 
Goldstein et al . 3 and Goldstein and Leib, 4 in that the viscous shear flow that 
bounds the nonlinear region is governed by the three-dimensional boundary 
layer equations. However, these equations now have to be solved on a curved 
surface, rather than on a plane surface, and this again makes the problem more 
complex. 

For definiteness, we suppose that the imposed spanwise sinusoidal 
crossflow has the form of a pair of counter-rotating vortices per period, 
which then become distorted as the flow evolves downstream. The numerical 
solution shows that the vortex sheet grows monotonically until the 
computations break down at (or near) ex = 0.8. The relevant Fourier spectrum 
suggests that this breakdown is due to the formation of a pair ot 
symmetrically disposed singularities per period. 

The breakdown is not completely unexpected, since the vortex sheet 
boundary condition has the form of an inviscid Burger's equation, with some 
additional terms arising from the vortex sheet development, and the Burger's 
equation solutions are known to become singular at finite downstream 
locations. The present problem also bears some similarity to the (temporally 
evolving) Kel vin-Helmhol tz instability of a disturbed two-dimensional vortex 
sheet, which is known to produce a finite time singularity. 5 Krasny 6 obtained 
a numerical solution to this problem and found that when the initial 
disturbance is sufficiently large, a pair of singularities ultimately form and 
the attendant Fourier spectrum is similar to the one obtained in the present 


3 


case. Krasny 6 used an ad hoc approximation to extend his computations beyond 
the singular time and found that the vortex sheet rolls up into a pair of 
tightly wound spirals. We therefore expect that the vortex sheet will also 
undergo spanwise roll-up in the present case. 


The shear layer flow appears to be relatively unaffected by the vortex- 
sheet distortion. The layer thickness does increase faster than the 

corresponding Lock profile, but there does not seem to be any dramatic changes ' 

I 

as the singularity is approached. This may reflect the fact that the j 

singularity is of higher order, i.e., it only appears in the higher-order 
derivatives. 


2. Formulation 

Consider a uniform, parallel, incompressible, steady stream that slips 
past a semi -infinite but infinitesimally thin flat plate with speed U„. A 
Cartesian coordinate system is centered at the plate's trailing edge, with the 
x-axis aligned with the main stream and y and z coordinates in the transverse 
and spanwise directions, respectively. In the large-Reynolds-number limit 
being considered herein, the undisturbed downstream flow consists of an 
inviscid parallel flow with a plcne vortex sheet in x > 0, y = 0 and a viscous 
layer that resolves the resulting velocity discontinuity. Lock 1,2 obtained a 
similarity solution to the viscous flow. We suppose that a steady crossflow 
disturbance of 0(e) (e « 1), with spanwise lengthscale X, is superposed on 
this piecewise uniform flow. All quantities are nondimensional ized using 
appropriate combinations of X and U,,. The resulting Reynolds number 
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R = XU«/v, where v is the fluid's kinematic viscosity, is assumed to be 
sufficiently large (R » e' 1 ) so that to 0(e z ) the outer flow satisfies the 
Euler equations: 

V-u=0 , 

u-Vu+Vp = 0 , (1) 

where u is the velocity, and p is the pressure. (The 0(R' A ) boundary layer 
displacement effects can be incorporated into u if necessary.) The vortex 
sheet, located at y = f(x,z), is a material surface and satisfies a kinematic 
requirement that it move with the fluid velocity, which means that for the 
steady flow being considered here, the component of velocity normal to the 
sheet must vanish: 

u-V(y-f) = 0; x> 0 , y = f (x, z) , -°°< z < °° . (2) 

Also, the pressure must be continuous across the sheet in order to ensure that 
it does not undergo infinite acceleration. 

Since the base flow is zero in the region below the sheet, the 
corresponding induced flow, if it exists, must be irrotational to 0(e 2 ) and 
therefore possesses a velocity potential, say However, $' turns out to be 
a constant, because it is an everywhere bounded solution to Laplace's 
equation, whose normal derivative vanishes on the (stationary) vortex sheet 
and plate boundaries. Hence there is no induced flow, the pressure is 
everywhere constant in y < f(x,z), and the second boundary condition become, 

p = 0; x> 0 , y= f{x, z) , -°°< z<°° . (3) 

The velocities and pressure must remain bounded at large distances from the 
trailing edge, and the upstream initial conditions require 

u - u 0 as x - (4) 
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where u 0 (y,z) = (0,v 0 ,w 0 ) is the imposed crossflow. Continuity of the vortex 
sheet requires that 

f= = 0; x = 0 , -°°< z<°° . (5) 

dx 

2.1. Breakdown of linear solution 


Near the trailing edge, the dependent variables expand like 

u(x,y, z) = (1,0 ,0) +e(u 1 ,v 1 ,w 1 ) +e 2 (u 2 ,v 2 ,w 2 ) +■■■ , (6) 

p(w,y,z) - ep 1 +e 2 p 2 +- , (7) 

f(x,z) = ef 1 +e 2 f 2 +- , (8) 


Substituting this into equations ( 1 ) - (3 ) and equating powers of e, we find 
that the 0(e) solution is f, = p n = 0 and 

{o, v 0/ w 0 ) , (y) 


while the 0(e 2 ) solution satisfies 

V'llj = 0 , 


dl 

dx 


+ Vp 2 = - (Uj/V) Ui 


(10) 

(ID 


subject to 


df 2 

dx 


= 0; x>0,y = 0, -<*><z< 00 , 


p 2 = 0; x> 0 ,y= 0, -oo<z<». 


(12) 


i 


i 

i 

i 

I 


p 2 must remain bounded in the far field, and taking the gradient of 
equation (11) and using equation (10) shows that p 2 satisfies the Poisson equation 
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(13) 


V 2 p 2 = -V* [ (Uj/V) uj , 

while the y-component of equation (11) implies that 

= 0; xzQ , y = 0 , -°°< z<°° . (14) 

dy 

Once p 2 is found from these relations, u 2 is a simple quadrature of (10). 

It follows that 

p 2 (x, y, z) ~p 2 „(y f z) {x~' a> ) , (15) 

since both the forcing in equation (13) and the boundary condition (14) become 
independent of x, as x -* Equations (10) and (11) imply that the 
corresponding asymptotic behavior of the velocity and sheet displacement are: 

u 2 — [ (Uj/V) Ui + Vp 2 J x (x-~) , (16) 

f 2 ~- (ivV) v i*~^ (^-*°°) • (17) 

This shows that 0(e 2 ) terms become comparable to the 0(e) terms when ex ~ 0(1) 
and, therefore, that the expansions (6)-(8) are no longer valid in this 
region. 

2.2 Nonlinear inviscid region 

To find an appropriate solution for this region (which will turn out to 
be uniformly valid when x = 0(1)), we introduce the slow streamwise length 
scale 

x = ex (18) 
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and seek an expansion of the form 

u = l+e 2 u(x,y, z , ) +••■ v=ev(x, y, z) +••• w=ew(x,y, z) +■•• p + e 2 p(x,y, z) 

f= fix, z) +••• . 


Substituting these into equations (1) - (3) and equating powers of e gives 


ix + 

dy 



( 20 ) 


Du 


_ dp 

ax' 


Dv = 


3p 

ay' 




9p 

az' 


( 21 ) 


where 


D 




( 22 ) 


and 


= 0; x> 0 ,y=0, -°°< z<°°, (23) 

ax az 


p = 0; x> 0,y=0, -°°< z<«>. (24) 

Matching with the perturbation solution (6) of §2.1 requires that 

v,w~v Q ,w 0 as x-0. (25) 


Eliminating p between equations (21) and (22) yields 

nQ= 0 (26) 


where 


o = OiL - — 

- dy dz 


(27) 


is the streamwise vorticity. The simple (Prandtl) coordinate transformation 

y = y- f (x, z) , (28) 
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transfers the vortex-sheet boundary to the known location y = o. The 
vorticity transport equation (26) remains 


but now with 


DQ = 0 


D = 


_1 + 

dx 



( 29 ) 


(30) 


where 


v = v-f--wf z . 


(31) 


The kinematic condition (23) is now simply 

v=0 at y = 0 . (32) 

Equation (20) implies that there exists a crossflow stream function iji 
defined by 




(33) 


Inserting this together with equation (31) into definition (27) shows that 

Q = U + f 2 z ) ^ w -^yf zz -2^ 7z f z -7- z . (34) 

Transforming the second two equations (21) via equations (28) and (31), 
eliminating p from the result, and using conditions (24) and (32), we find 
that the latter two conditions imply 

D 0 w+I z D*f=0; at y = 0 , (35) 

where 
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Thus the problem consists of finding a solution to equations (25), (29), and 
(33) through (35) that remains bounded at y-°°. 


2.3 Viscous layer 

The vortex sheet discontinuity is resolved by a thin viscous shear layer 
of thickness, 5 = (L/X)0(R l ' % ) , where R L is the Reynolds number based on the 
length L from the origin of the boundary layer, so that 6 = (eR)' % when 
x = 0(1). We therefore introduce the scaled transverse coordinate, Y, defined 
by 

y = bY. (36) 


Then the streamwise velocity 

u — U(x, Y, z) (37) 

is 0(1), and the requirements of continuity and matching with the outer flow 
are ensured by the scalings, 

e6v = v-euf^-wf. t , w = eW, p = e 2 P. (38) 

Introducing these new scaled variables into the Navier-Stokes equations and 
retaining only lowest order terms, we find 

DU = (1 + fg)U YY , (39) 

( 1 + f]) DW+ f z ( UDf-z + WDf z ) = (1 + f\) 2 ^yy , (40) 

Ujt+V Y +W z = O , (41) 


where 


0ma i 


W- 


dz 


The boundary conditions are obtained by matching to the outer flow, 
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U-* 1 , W-W m s w(.x,y = 0 , z) (r-«) 

U- 0 , W- 0 (y--«) , 


( 42 ) 


Equation (40) becomes identical to (35) when all derivatives with respect to Y 
are set equal to zero, as required for smooth matching with the outer flow. 

The initial condition is Lock's 1 solution which is valid for x = 0(1), i.e., 
for x - 0 . 


3. Results 


In order to be specific, we consider only the imposed crossflow 


v= coz z sech y tanh y , 
iv=-sin z sech y{2 tanh 2 y- 1) , 

Q = -6 sin z sech y tanh y(tanh 2 y-l) . 


which is initially sinusoidal and decays exponentially as y-°°. Both the 
inviscid and viscous flow can be computed by marching downstream in x. 

3.1. Inviscid flow 

Since the initial crossflow was restricted to spanwise periodic data, 
Fourier spectral methods can be used in spanwise direction, and we use second- 
order, finite-difference formulae on uniform grids to compute streamwise and 
transverse derivatives. The vorticity is obtained at each station x 1 by the 
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discretization of equation (29) according to the mid-point, leap-frog method, 
which is explicit and marginally stable. Next, the stream function 
i|r (x i ,y,z) and vortex-sheet displacement f(x 1 ,z ) are determined 

sequentially by iteration between equations (34) and (35). All nonlinear 
terms are evaluated as products in physical space. The computed solutions 
were tested against a perturbation solution for small x. The numerical 
results (presented below) were obtained using A5c=0 . 005, Ay=0 . 005 , with 
the computational region limited to y<. 5. . 

The vortex sheet in the region 0 <x< 0 . 8, 0 < z< 2n is shown in figure 2. 
Its spanwise form is dominated by the second harmonic of the imposed 
crossflow. Figures 3 and 4 show the downstream development of sheet 
displacement f(x,z) and the spanwise velocity w(x,y=o,z) at several 
intermediate stations. The steepening of the velocity profile is reminiscent 
of the solutions to the inviscid Burger's equation which equation (35) 
resembles. However, this equation also contains additional terms, due to 
vortex-sheet displacement which reverses the steepening (about z = n) near 
x=0.l . This can be seen more clearly in figure 5, which is a plot of the 
derivative, dw/dz. Once the steepening reverses, two symmetrically disposed 
lobes emerge and continue to sharpen until the computations break down. 

Figure 6 is a contour plot of a crossflow steam function if' defined 
analogously to i|i, but based on the Cartesian variables v,w,y,z. The 
vortices are attracted towards the z = n symmetry plane, and an induced 
rotation opposite to that in the main flow appears close to the vortex sheet. 
(Contours end on the vortex sheet because a curved sheet implies a non-zero 


12 



velocity in the crossflow plane.) These streamline patterns do not suggest 
any obvious tendency for singularity formation or breakdown. 

The breakdown occurs near 5c= 0.8. Calculations were repeated with the 
number of spanwise nodes doubled successively from 32 through 512, 
corresponding to increasing the highest resolved wavenumber from 16 to 256. 

In all cases, the breakdown occurs at some x=x g >0.8 with x s approaching 0.8 

from above, monotonically, as the resolution improves. Moreover, the 
breakdown occurs on the sheet and not at an interior point. 

The route to singularity formation can be observed in the downstream 
evolution of the spectrum of the solution. Figure 7 is the Fourier spe'.crum, 
Q n (x,y= 0) , of the spanwise velocity profiles shown in figure 4. The 

spectrum must exhibit exponential decay at large wavenumbers as long as the 
function remains analytic, but energy transfers to higher wavenumbers, and the 
decay rate decreases as x increases. Near x = 0.8, the exponential decay is 
eventually lost, and the function is no longer analytic. The nearly periodic 
oscillations in the high-wave-number component of the spectrum suggest that 
the singularity initially forms at two spanwise locations, say z = ±<. s per 
period. The interval An between successive minima in the spectrum is 
approximately 16, suggesting that z s = tt/ 16, which is approximately the 
distance from the symmetry plane of the two minima in figure 5. The spectrum 
does not seem to suggest that the oscillatory tail is due to spurious growth 
at the highest wavenumbers. 
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3.2. Viscous layer 


The viscous layer is driven by the tangential velocity over a distorted 
vortex sheet. It is computed by an adaption of the Keller-box scheme, 7 which 
is used to march downstream from the Lock profile. 2 

The layer thickness is expected to grow, and the governing equations 
(42) - (44) are rewritten in terms of the Blasius similarity variable q=Y//x. 

The box scheme allows arbitrary spacing in both x and rj and is second-order 
accurate since the discretized equations are applied at the centers of cells 
formed by the x-q grid. The nodes were distributed so that the spacing 
between nodes increases at a fixed rate towards the upper and lower edges «f 
the viscous layer: 

V = /(l -h 2 ) (j = 1,2, ...<7*) ; 

here h 1 = 0.005 and h 2 = 1.05. The range -21 < q < 13 corresponding to 
J* = 110 and J + = 100 was found to be adequate for the inner solution to 
smoothly reach its asymptotic boundary values. While spectral techniques were 
used conveniently in the inviscid region, they become inordinately expensive 
as the resolved spanwise wavenumber is increased beyond 32 (many general 
matrices of the order of the maximum wavenumber need to be inverted), and 
instead, second-order, l inite-difference, upwind formulae were used in the 
viscous solution. The symmetry of the external flow is imposed on the viscous 
solution and, as is well known, the solution on the planes of symmetry, 
z = 0 ,tt , can be obtained independently of the solution in the region between 
these planes. We therefore calculated the solution on the z = 0 symmetry 
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plane and obtained the solution on interior lines by marching from z = 0 to 
z = jt. Since the spanwise velocity W remains positive everywhere in 
0 < z < 7r, this can be carried out without iteration. 

The external crossflow and its distortion of the vortex sheet cause the 
flow to depart from the Lock profile, with the most noticeable changes being 
the increase in layer thickness on the z = it plane and the decrease on the 
z = 0 plane. Figure 8 shows the increasing distortion of the streamwise 
velocity profiles with increasing downstream distance. Figure 9 shows the 
corresponding spanwise velocity profiles. Vector plots of the crossflow are 
shown in figure 10. Note that the transverse velocity is 0(eR)‘ 1/2 smaller 
than the spanwise velocity. The pair of transverse jets (figure 10) at the 
upper edge of the viscous layer correspond to the reversal in the steepening 
of the spanwise velocity in figures 4 and 5. This follows from expanding the 
inviscid solution in a Taylor series and using continuity to obtain 



4. Concluding Remarks 

Plane mixing-layer flows, being inviscidly unstable, have been observed 
to rapidly break down and roll up into tightly wound spirals. Two-dimensional 
mixing effectiveness is inadequate for many technological applications, and 
three-dimensional mechanisms, such as spanwise roll-up, are sought to increase 
mixing. In this paper, we show that the cumulative action of weak spanwise 
vorticity in the external free stream can lead to an 0(1) distortion of the 


15 


mixing region. Since the vorticity, in our analysis, resides in the external 
flow, it can be maintained over large streamwise distances without undergoing 
viscous decay. The formation of an inviscid singularity in the solution 
suggests, by analogy with studies of the similar Kelvin-Helmholtz problem, 
that spanwise roll-up will ultimately occur in this flow. 

Even though the analysis cannot be extended beyond the initial 
singularity, it can still yield some useful information about mixing by using 
it to investigate the effect of the initial vorticity distribution on the 
location of singularity formation. Early formation of the singularity should 
lead to enhanced overall mixing. 

The authors would like to thank Pr. Stephen J. Cowley of the Department 
of Applied Mathematics and Theoretical Physics, University of Cambridge, for 
his helpful suggestions about the singularity and Dr. Stewart Leib for 
providing us with his computer program. 
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Figure 1 - Flow Configuration 
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Figure 8 (a) Strear.wse Velocity x = 0.3 
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Figure 8 (b) Stream. wise Velocity 
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Figure 8 (c) Strearrvvise Velocity x = 0.795 
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Figure 9 (a) Spa-.wise Velocity x = 0.3 
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Figure 9 (b) Spanwise Velocity x = 0.7 
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Figure 10 (a) Crossflow Velocity Vector 
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